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Over the last two decades, both experiments and simulations have demonstrated that 
transverse wall oscillations with properly selected amplitude and frequency can reduce 
turbulent drag by as much as 40%. In this paper, we develop a model-based approach 
for designing oscillations that suppress turbulence in a channel flow. We utilize eddy- 
viscosity-enhanced linearization of the turbulent flow with control in conjunction with 
turbulence modeling to determine skin- friction drag in a simulation- free manner. The 
Boussinesq eddy viscosity hypothesis is used to quantify the effect of fluctuations on the 
mean velocity in the flow subject to control. In contrast to the traditional approach that 
relies on numerical simulations, we determine the turbulent viscosity from the second 
order statistics of the linearized model driven by white-in-time stochastic forcing. The 
spatial power spectrum of the forcing is selected to ensure that the linearized model for 
the uncontrolled flow reproduces the turbulent energy spectrum. The resulting correction 
to the turbulent mean velocity induced by small amplitude wall movements is then used 
to identify the optimal frequency of drag reducing oscillations. In addition, the control 
net efficiency and the turbulent flow structures that we obtain agree well with the results 
of numerical simulations and experiments. This demonstrates the predictive power of our 
model-based approach to controlling turbulent flows and is expected to pave the way for 
successful flow control at higher Reynolds numbers than currently possible. 



1. Introduction 

1.1. Background 

Turbulent flows are ubiquitous in nature and engineering. Dissipation of kinetic energy 
by turbulent flow around airplanes, ships, and submarines increases resistance to their 
motion. This motivates design of control strategies for enhancing performance of vehicles 
and other systems involving turbulent flows. Utility of different approaches for maintain- 
ing the laminar flow, reducing skin-friction drag, and preventing separation is surveyed 
in Joslin (1998); Gad-el Hak (2000). While traditional flow control techniques combine 



physical intuition with costly numerical simulations and experiments, model-based tech- 
niques utilize developments from control theory to improve flow manipulation. Recent 
research suggests that traditional strategies can be significantly enhanced by flow control 
design based on analytical models and optimization tools ( Kim fc Bewley||2007 ). 

The effectiveness of model-based feed back (Hogberg, Bewley fc Henningson|2QQ3j ) and 
sensor-less ( Moarref fc Jovanovic||2010 Lieu, Moarref fc Jovanovic||2010 ) techniques for 
controlling the onset of turbulence at low Reynolds numbers stems from their ability 
to reduce receptivity and enhance robustness of the flow. Model-based approach to flow 



2 



R. Moarref & M. R. Jovanovic 



control design has been motivated by realization that a mechanism which initiates tran- 
sition is governed by the degradation of rob ustness (Trefethen et a/. [1993 |Schmid 2QQ7| 



and the consequential noise amplification (Farrell & loannou 1993[ Bamieh & Dahleh 



2001 Jovanovic 2004 Jovanovic & Bamieh"'2005). Consequently, the above mentioned 



techniques have utilized Navier-Stokes (NS) equations linearized around the laminar flow 
as a control-oriented model with the objective of reducing sensitivity to modeling imper- 
fections. In addition, a fully-developed turbulent channel flow at low Reynolds numbers 
was relaminarized via a gain-scheduled linear state-feedback controller (Hogberg, Bew^ 
ley fc Henningson|2003 ). Comparison between the impulse responses of the NS equations 



linearized around the laminar and turbulent mean velocities and the direct numerical sim- 
ulations (DNS) of the turbulent flow subject to small amplitude impulsive perturbations 
has been provided by Luchini, Quadrio & Zuccher (2006). 



The role of linear mechanisms in formation and maintenance of streamwise streaks in 
turbulent flows was examined by Lee, Kim fc Moin^ (1990 J . It was shown that the streaks 



are formed by linear amplification of eddies that interact with the large mean shear. 
For homogenous flows subject to high shear rates, Lee et al. (1990) demonstrated that 
the DNS-based instantaneous velocity is similar to the velocity that is predicted by the 
linearized equations in the rapid distortion limit ( Pope|[2QQQ ). Furthermore, Kim fc Lim 



(2000) used DNS to demonstrate the importance of the linear vortex tilting mechanism 



in maintaining the streamwise vortices in a fully-developed turbulent channel flows. 



Schoppa fc Hussain (2002p conducted a thorough investigation of the role of near- wall 



streaks in generation of streamwise vortices. By examining the evolution of infinitesi- 
mal fluctuations around a streaky turbulent base flow, they determined the amplitude of 
streaks above which modal instability occurs. It was noted that only 20% of the streaks 
in the fully developed turbulent buffer layer are strong enough to trigger modal insta- 
bility. Consequently, a secondary transient growth mechanism for bypass transition in 



wall-bounded shear flows was proposed (j Schoppa fc Hussaiii||2002 Hoepffner, Brandt fc 
|Hen ningson"2005') ; this mechanism appears to be capable of producing much larger tran- 
sient growth rates than the secondary streak-instabilities. Furthermore, , Chernyshenko fc] 



Baig ( 2005| ) studied the origin of near- wall streaky patterns and demonstrated that the 



combination of linear effects including lift-up of the mean profile, tilting and stretching 
by the mean shear, and viscous dissipation induce formation of these patterns. Capability 
of the linearized NS equations to qualitatively predict both the streak spacing and its 
dependence on the wall-normal distance was also shown. 

1.2. Previous studies on drag reduction by transverse wall oscillations 

Several experimental and numerical studies have shown the effectiveness of sensor-less 
strategies for turbulence suppression in wall-bounded shear flows. The experiments of Brad- 
(1985) and DNS of Moin et al. (1990) showed that imposing a constant 



shaw fc Pontikos 



transverse strain on a turbulent boundary-layer can transiently reduce the turbulent 



kinetic energy and the Reynolds stresses. Motivated by this observation, Jung, Man- 



giavacchi fc Akhavan (1992) used DNS to establish a sustained turbulence suppression in 
a channel flow subject to transverse wall oscillations. For the flow with friction Reynolds 
number = 200, skin-friction drag reduction of up to 40% was reported with maxi- 
mum drag reduction taking place for the period of oscillations T+ ^ 100 (in viscous time 



units). The numerical results of Jung et al. 


(1992) were experimentally verified by|Laad- 


hari, Skandaji fc Morel ( 


1994); 


Choi, DeBisschop fc Clayton 


(1998); ( 


C!hoi| (|2002 


) ; |Ricco 


(2004), where a drag reduction of up to 45% was reported. Choi 


(200^ 


I ) argued that wall 



oscillations induce negative spanwise vorticity, thereby suppressing turbulence by ham- 



pering the vortex stretching mechanism. In addition, the experiments of Ricco (2004) 
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showed that the near-wall flow is dragged laterally by wall oscillations which reduces 
the length of the streaks and increases the spacing between them. [Touber fc Leschziner| 
(2012) also demonstrated that wall oscillations significantly distort the near- wall streaks 
and reduce the contribution of turbulence to the wall shear stress. Recent DNS study 
of Ricco et al. (2012) further revealed that wall oscillations directly affect the turbulent 
dissipation. For T+ < 100, it was shown that drag reduction scales linearly with the 
volume integral of an enstrophy production term caused by the spanwise shear layer. 

Several alternative mechanisms for inducing transverse oscillations have also been in- 
vestigated. For example, Berger et al. (2000) used DNS of conductive fluids in a channel 
flow with Rr = 100 to show that time-periodic spanwise Lorentz force can reduce skin- 
friction drag up to 40%. The amount of drag reduction was found to decrease for larger 
Rr. In addition, Du & Karniadakis (2000); Du, Symeonidis & Karniadakis (2002) studied 
the effect of Lorentz force in the form of spanwise traveling waves confined to the viscous 
sub-layer. For R^ = 150, their DNS showed a drag reduction of up to 30%. The drag- 
reducing mechanisms of transverse motions induced by spanwise oscillations, spanwise 
traveling waves, and riblets have been surveyed by Karniadakis & Choi (2003). Recently, 
turbulent drag reduction by waves of spanwise velocity that travel in the streamwise direc- 
tion has been examined using DNS ( Quadrio, Ricco fc Viotti|2009 ), experiments (Auteri 



et a/.|[2QTQ ), and generalized Stokes layer theory ( Quadrio fc Ricco||2011 ). It was shown 



that upstream traveling waves reduce drag at any speed. On the other hand, down- 
stream waves reduce drag only at speeds that are much larger or much smaller than the 
convecting speed of near- wall turbulent structures. 



Quadrio fc Sibilla (2000) used DNS to show that up to 40% of drag reduction can 



be achieved by oscillating a cylindrical pipe along its longitudinal axis. In a series of 
papers, [Baron fc Quadriol ( |1996D ; 'Quadrio fc Ricco] ( |2003[ [20041 ); [Ricco fc Quadrio] ( [2008 ) 
further studied drag reduction by transverse wall oscillations in a channel flow; also see 
recent review by Quadrio (2011). In addition to quantifying the saved power associated 
with drag reduction, they accounted for the input power necessary for maintaining wall 
oscillations; for small oscillation amplitudes, it was established that a net power gain with 
drag reduction of up to 10% can be achieved. Furthermore, for the same values of oscilla- 
tion amplitude and frequency, DNS of Choi, Xu fc Sung (2002) showed that the amount 
of drag reduction in a channel flow can drop by as much as 25% as the friction Reynolds 
number R^ increases from 100 to 400. More recent DNS studies of [Ricco fc Quadno] 
(2008) and Touber fc Leschziner (2012) confirmed deterioration in drag reduction by 
15% and 18% with increase in R^ from 200 to 400 and from 200 to 500, respectively. 
On the other hand, for bounda ry layers subject to plate oscillations with small periods 
(T+ < 83), the experiments of Ricco fc Wu (2004) showed a weak dependence of drag 
reduction on the Reynolds number. Ricco fc Quadrio[ (^2008 ) argued that this discrepancy 
may arise from high experimental uncertainty, effects of geometry, and smaller values of 
oscillation periods relative to DNS studies. 

In contrast to the aforementioned experimental and numerical studies, there are rela- 
tively few theoretical developments regarding drag reduction by wall oscillations. We next 



briefly summarize notable exceptions (Dhanak fc Si 1999 Bandyopadhyay 2006 Ricco 



fc Quadrio|2008 ). Dhanak fc Si (1999) used exact solutions of the NS equations to study 



the interactions between the evolving streamwise vortices in a turbulent boundary layer 
and the Stokes layer induced by wall oscillations. Simulations of the resulting dynam- 
ical model showed that these interactions reduce the Reynolds stresses, the turbulence 
production, and the skin friction drag. Bandyopadhyay (2006) proposed a vorticity reori- 
entation hypothesis and showed that the negative spanwise vorticity induced by the wall 
oscillations modifies the orientation of the total vorticity field in the near wall region 
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and suppresses turbulence production. Furthermore, it was shown that the developed 
model yields drag reduction that agrees well with experimental data. Finally, in addition 



to conducting a thorough DNS study, Ricco & Quadrio (2008) used the solution to the 



laminar Stokes problem to quantify the dependence of the drag and the input power 
necessary for maintaining wall oscillations on their period and amplitude. 

1.3. Preview of key results 

While most model-based efforts to date have considered the problem of maintaining the 
laminar flow and relaminarization, in this paper we show that turbulence modeling in 
conjunction with eddy- viscosity-enhanced linearization can extend utility of these meth- 
ods to control of turbulent flows. Control-oriented turbulence modeling is challenging 
because of complex flow physics that arises from strong interactions between the turbu- 
lent fluctuations and the mean velocity. We build on recent research that demonstrates 
considerable predictive power of nonmodal stability analysis (Schmid__2007) even in tur- 



bulent flows (del Alamo & Jimenez 



2006 



Cossu et al. 



2009 



Pujals et a/.|2009). These 



papers have shown that the equations linearized around turbulent mean velocity, with 
molecular viscosity augmented by turbulent viscosity, qualitatively capture features of 
turbulent flows with no control. For the flow with control, we examine the class of lin- 
earized models considered by 



del Alamo & Jimenez 



(2006); 



Cossu et al. 



(2009); Pujals 



et al. |2009 ) and use turbulent viscosity hypothesis to quantify the influence of turbu- 



lent fluctuations on the mean velocity. We demonstrate the ability of this approach to 
quantitatively predict the effect of control on turbulent drag. 

The difficulty here arises from the fact that the turbulent viscosity of the flow with 
control has to be determined. Even though we use the Boussinesq hypothesis to capture 
the influence of control on turbulent viscosity, in contrast to current practice we do not 
rely on numerical simulations for flnding turbulent kinetic energy k and its rate of dissi- 
pation e. Instead, we introduce a simulation-free method based on stochastically-forced 
linearized model of controlled flow to obtain k and e from the second-order statistics of 
velocity fluctuations. These statistics are used to determine the turbulent viscosity for 
the flow with control, and thereby to compute the effect of control on the turbulent mean 
velocity and on the skin-friction drag. 

We utilize linearized equations subject to white-in-time stochastic forcing with appro- 
priately selected second-order spatial statistics. Using analogy with homogenous isotropic 
turbulence ( Jovano vic fc Georgiou|[2QTQ ), we select these to be proportional to the two- 
dimensional energy spectrum of the uncontrolled flow. Note that while our approach 
takes advantage of the turbulent viscosity and the energy spectrum resulting from direct 
numerical simulations { DNS) of the uncontrolled flow ( Kim, Moin fc Moser|1987 _ Moser, 



Kim & Mansour 



1999 



del Alamo & Jimenez 



2003 



del Alamo, Jimenez, Zandonade fc 



Moser|2QQ4 ), we do not rely on numerical or experimental data for determining the effect 



of control on the turbulent flow. To the best of our knowledge, the present work is the flrst 
to utilize publicly available DNS data of the uncontrolled flow to guide control-oriented 
modeling of turbulent flows. Even though the aforementioned databases provide coarse 
one-point correlations in the wall-normal direction, we demonstrate that they can be 
effectively employed for model-based flow control design. 

In this paper, we use a model-based approach to examine the effect of transverse wall 
oscillations on the dynamics of a turbulent channel flow. We start by showing that the 
power necessary for maintaining wall oscillations increases quadratically with their am- 



plitude, which is in agreement with DNS of Quadrio fc Ricco ( 2004p and theoretical study 
of Ricco fc Quadrio (2008). Since large control amplitudes yield poor net efficiency, we 
confine our study to small oscillation amplitudes and use perturbation analysis (in the 
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Figure 1. (a) Pressure driven channel flow; and (b) Channel flow subject to transverse wall 

oscillations. 

amplitude of oscillations) to identify the period of oscillations that achieves largest drag 
reduction in a computationally efficient manner. In addition, we quantify the net effi- 
ciency, discuss the drag reduction mechanisms and the effects of the Reynolds number, 
and compare the dominant structures in flows with and without control. The close agree- 



ment between our results and the results obtained in experiments and DNS (Jung et al. 



1992| [Baron fc Quadrio||1996[ |Choi||2QQ2| [Quadrio fc Ricco||2QQ4[ [Ricco fc Quadrio||2QQ8 ) 

demonstrates the predictive power of our model-based approach to flow control design. 

Our presentation is organized as follows: in § [2j we formulate the problem and provide 
a brief overview of the governing equations, turbulent mean velocity, skin- friction drag 
coefficient, net efficiency, and the Boussinesq eddy viscosity hypothesis. In § |3j we use 
eddy-viscosity-enhanced stochastically forced linearized model to study the dynamics of 
infinitesimal fluctuations around the turbulent base flow. We also describe an efficient 
method for computing the second-order statistics of fluctuations. These statistics are 
used to determine the influence of control on turbulent viscosity. In § [4j we apply our 
theoretical developments to the problem of turbulent drag reduction with transverse wall 
oscillations, and provide a thorough analysis of the effect of control on skin-friction drag 
and net efficiency. In § [5) we use characteristic eddy decomposition to visualize the effect 
of control on turbulent flow structures. We conclude the paper with a brief summary of 
our contributions and outlook for future research in § [6j 



2. Problem formulation 

The pressure-driven channel flow of incompressible Newtonian fluids, with geometry 



shown in figure 1(a), is governed by the non-dimensional NS and continuity equations 

u, = -(u • V)u - VP + {l/Rr) Au, 
= V-u, ^ ' ^ 

where u is the velocity, P is the pressure, V is the gradient, and A = V • V is the 
Laplacian. The Reynolds number = Urh/u is defined in terms of the channel's half- 
height h and the friction velocity u^- = ^Jr^^j (x, z) are the streamwise, wall-normal, 
and spanwise coordinates, and t is time. Here, is the wall-shear stress (averaged over 



horizontal directions and time), p is fluid density, and v is kinematic viscosity. In (2.1) 
and throughout the paper, spatial coordinates are non-dimensionalized by /i, velocity 
by time by h/ur^ and pressure by pu^. When normalized by v/ur^ the wall- normal 
coordinate is denoted = Rr{l ^ y). The subscripts x, y^ z, and t are used to denote 
the spatial and temporal derivatives, e.g., Ut = du/dt = dtu. 

Throughout the paper we assume that the bulk flux, which is obtained by integrating 
the streamwise velocity over spatial coordinates and time, remains constant. This con- 
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straint is commonly imposed in experiments and DNS of turbulent flows and it can be 
enforced by adjusting the uniform streamwise pressure gradient P^. In addition to the 
driving pressure gradient, which balances the wall-shear stress ( McComb||1991[ ), the flow 
is also subject to zero- mean transverse wall oscillations of amplitude a and frequency cjt; 
see figure [T(b)] The period of oscillations normalized by h/ur (outer units) is given by 
T = 27r/uJt, which is equivalent to T+ = R^T when normalized by i^/u'^ (viscous units). 
The streamwise and wall-normal velocities satisfy no-slip and no-penetration boundary 
conditions at the walls. 

Reynolds decomposition separates the velocity in a turbulent flow into the sum of the 
turbulent mean velocity, \J = [U V W]^^ and the fluctuations around U, 



V = \ U V w \ 



u = U + V, U = (u) , (v) = 0. 
Here, ( • ) denotes the expectation operator, 

1 f'f 



{u{x,y,z,t)) 



1 

lim — / u(x, t + r) dr. 

tf^ootf Jo 



This decomposition yields the Reynolds-averaged Navier-Stokes (RANS) equations for 
the turbulent mean velocity (McComb„1991, Durbin fc Reif 2000} |Pope||2QQQD , 

U, = - (U.V)U - V(P) + (l/i?.)AU - V.(vv^>, 

= v-u. 

In a turbulent flow, the second-order statistics of fluctuations (vv-^), i.e., the Reynolds 
stresses, introduce additional flux. The Reynolds stress tensor quantifies the transport of 
momentum arising from turbulent fluctuations and it has profound influence on the mean 
velocity, and thereby on skin- friction drag ( ^McComb|1991 ). The difficulty in determining 
statistics of fluctuations arises from the nonlinearity in the NS equations which makes 
the nth velocity moment depend on the (n + l)th moment ( McComb||1991 ). 



2.1. The turbulent mean velocity 



The closure problem in (2.2) can be overcome by expressing the higher order moments 



in terms of the lower-order moments. According to the Boussinesq eddy viscosity hy- 
pothesis ( |McComb|ri991[ [Durbin fc Reif|[2QQQ{ |Pope||2QQQD , the turbulent momentum is 
transported in the direction of the mean rate of strain. 



itr((vv-} 



^ (VU + (VU)^) 



(2.3) 



where i^riy) is the turbulent viscosity normalized by z^, overline denotes averaging over x 
and tr (•) is the trace of a given tensor, and / is the identity tensor. In the flow subject 
to wall oscillations, ut is a function of both y and t. However, in this paper, we consider 
the averaged effect of control (over one period of oscillations T) on ut- This does not 
mean that the time-periodic flow is treated as time-independent. In fact, the approach 
outlined in § |3.1| and in Appendix [C] facilitates analysis of the time-periodic flow. 

The steady-state solution (after the influence of initial conditions disappears) of the 
system (2.2)- (2.3) subject to a uniform pressure gradient, P^^, and the transverse wall 



oscillations, 

W{y ±l,t) = 2asm{ujtt) , 

is determined by [U{y) W{y^ ^) It can be shown that the streamwise mean velocity 
averaged over one period of oscillations, U{y), arises from the uniform pressure gradient. 
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and that the spanwise mean velocity, W{y,t)^ is induced by the wah oscihations, 

r = {{i + My))Uy{y))y - RrP., 

{ U{y = ±1) = 0, 

r RrWt{y,t) = {{l + My))Wy{y,t))^ 

\ W{y = ±l,t) = 2asm{L0tt). 

Here, (1 + vt) represents an effective viscosity that accounts for botfi molecular and 
turbulent dissipation ( Pope|pOOO ) . 



(2.4a) 
(2.46) 



For given turbulent viscosity vt and driving pressure gradient P^-, \2A\ is an uncoupled 

ed 



system of equations for U and W; U can be obtained by solving (2.4a) 

U{y) = RrPx ^ , " ,., dC, 



(2.5) 



and W can be obtained by solving (2.46). The imposed boundary conditions induce a 



^/-dependent time-periodic spanwise velocity, 

W{y,t) = a(M^p(y)e"^'* + W^;(y) e"^'^* *), (2.6) 

where i = >/— 1, and * denotes the complex conjugate. An equation for Wp{y) is obtained 
by substituting ( p^ in ( |2.4&| ) 

\RrUJtWp{y) = {i + My))w^{y) + My)K^y), 

%(±1) = -i, 



(2.7) 



where prime represents differentiation with respect to y. The expression for Wp{y) is 
readily obtained from a solution to the resulting two point boundary value problem (2.7). 
The difficulty in determining U and W from (2.4) arises from the fact that vt depends 



on the fluctuations around the turbulent mean velocity, and thus it is not known a priori. 
A significant body of work has been devoted to finding an expression for tjt that yields the 
turbulent mean velocity in the uncontrolled flow (Malkus| |1956 Cess"1958 Reynolds & 
Tiederman|1967 ) . Reynolds fc Tiederman (1967) extended the turbulent viscosity model, 
originally introduced by Cess (1958), from the pipe flow to the channel flow. 




-il-\y\)RJci 



1/2 



(2.8) 



This expression employs the law of the wall in conjunction with van Driest 's damping 



function and Reichardt's middle law (Reynolds & Tiederman 1967). The parameters 



ci and C2 appear in the van Driest 's wall law and in the von Karman's log law (Pope 
2000), respectively. These two parameters are selected to minimize least squares deviation 
between the mean streamwise velocity (|2.5|) obtained with = —1 and turbulent viscos- 



ity (2.8), and the mean streamwise velocity obtained in experiments or in simulations. 
Application of this procedure yields {Rr = 186, ci = 46.2, C2 = 0.61; i^^- = 547, ci = 
29.4, C2 = 0.45; Rr = 934, ci = 27, C2 = 0.43} for the corr esponding mean velocities in a 



turbulent channel flow resulting from DNS (del Alamo fc Jimenez|2003 



f- 
J2C 



del Alamo et al. 



2004D . 

Under the assumption that the turbulent viscosity (2.8) captures the effect of back- 



ground turbulence on the mean velocity, the system of equations (2.4)-(2.8) yields a 
solution Uo = [^o(^) Wo{y^t)]^ . By construction, Uq approximates the mean stream- 
wise velocity in the uncontrolled turbulent flow, and Wq is the spanwise velocity induced 
by the wall oscillations and obtained under the assumption that the turbulent viscosity 
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Figure 2. (Color online available atjournals.cannbridge.org/flm) (a) The real part; and (b) the 
imaginary part of the solution Wp^o{y~^ ;T~^) to ( |2.7| (solid curves) in the flows with Rr = 186 
and T+ = 30 (blue), 100 (green), and 300 (redjTThe solution Wp,o{y^;T^) is obtained un- 
der the assumption that the turbulent viscosity of the uncontrolled flow captures the effect of 
fluctuations on the mean velocity. The solutions of the laminar Stokes problem are shown for 
comparison (dashed curves). 



is not modified by control. As discussed below, this strong assumption is only used as a 
starting point for our analysis. Figure [2]shows that Wp^o{y^;T^) in flow with Rr = 186 
(solid curves) is localized in the viscous wall region, y^ < 50. As expected from the anal- 
ogy to the solution of the Stokes second problem (for example, see Panton ( 1996| )), Wp^o 
shifts away from the wall as T+ increases. In addition, for small values of T+, we ob- 
serve close correspondence between Wp^o and the solution to the laminar Stokes problem 



(dashed curves) determined by Ricco & Quadrio (2008); both of these are in agreement 



with t he DNS results of [Quadrio fc Sibilla| ( |2QQQ| ); |Choi etal] (|2002|; [Quadrio fc Ricco 
(2003). For larger values of T+ , the discrepancy between these solutions increases. Ricco 
& Quadrio ( 2008 ) argued that in this case the rapidly- varying time evolution of the near- 
wall turbulent structures is not properly captured by the slowly-varying Stokes layer. 

Therefore, if the turbulent viscosity of the uncontrolled flow uto is used to model 
the oscillations induce Wq but have no impact on Uq (which, in this case, arises only from 
the uniform pressure gradient). The implications of this assumption for determining the 
skin-friction drag coefficient and the control net efficiency are discussed in § 2.2 where 



we demonstrate the necessity of accounting for the effect of control on the turbulent 
viscosity. 

2.2. Skin-friction drag coefficient and net efficiency 

As mentioned in § |2j the pressure gradient Px has to be adjusted in order to maintain 
the constant bulk flux, 

1 

= 2 y U{y)dy = const. 

Since t he skin- friction drag coefficien t is proportional to \Px\ and inversely proportional 
to ( [McComb[[T99T| [Panton[[l996[ ), 

Cf = 2\Px\/Ul 

for the flow with constant Ub, reduction (increase) in \Px\ induces drag reduction (in- 
crease). The change in the skin-friction coefficient relative to the uncontrolled flow is 
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100 200 300 



50 100 150 200 250 300 

r+ 

(a) 

Figure 3. (Color online) (a) (Solid curve) The required power, %nreq,o(r^), as a function of 
the period of oscillations T+ fo r the flow with Rj 18 6: (Symbols) DNS-based %nreq(T+;a) 




normalized by a at Rr = 200 jQuadrio & Ricco 20041 for oscillation amplitudes a = 2.25, 



a = 6, □; and a = 9, V; (b) The solid curve and the symbols show the dependence on l/vT+ 
for the same quantities as in figure 3(a) The red d otted curve shows the re quired power obtained 
from the solution to the laminar Stokes problem ( [Ricco Quadrio|2008 ) . 



given by 

%Cf = IQQ^i^^^-^ = 100(1 + J, (2.9) 

where the subscripts u and c denote the quantities in the uncontrolled and controlled 
flows, respectively. Thus, the control leads to drag reduction when P^ c > — 1- 

The drag reductio n induces sav ing in power (per unit area of the channel surface and 
normalized by pu"^) ( Currie|[2QQ3 ) 



Hsave = 2/7b(1 + 



x,c J 



Compared to the power required to drive the uncontrolled flow, = 2Ub^ the saved 
power is determined by the relative change in the skin- friction coefficient, 

%nsave = 100^ = 100(1 + P,,e) = %Cf. 

On the other hand, an input power is required for balancing the spanwise shear stresses 
at the walls ( Currie|[2003 l) . The required power exerted by wall oscillations expressed in 
fraction of the power necessary to drive the uncontrolled flow is given by (see Appendix|A]) 

2 

%nreq = 100-^Im(H^;|^^_^ - ^^Ui)' (2.10) 

where Im(-) denotes the imaginary part of a complex number. The net efficiency of 
control is quantified by the difference of the saved and required powers 

%nnet = %nsave " %nreq. 

Since the net efficiency is obtained from U and determining the turbulent mean 
velocities is essential for assessing the efficiency of wall oscillations. 

For the spanwise mean velocity Wq determined in § |2.1[ the required power grows 
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quadratically with a, 



req 



req,0 



100 



Figure 3(a) shows that %nreq,o(^^) decreases monotonicany with T+ and that smah 



values of require prohibitively large input power. The inset demonstrates close corre- 
spondence between %nreq,o (solid curve) and the DNS results of Quadrio & Ricco| ( 2004 ) 



(symbols). Using the solution of the laminar Stokes problem, Ricco fc Quadriof (2008) 
showed that the requi red po wer scales as a^/VT+. For T+ < 150 (or, equivalently, for 
l/v^T+ < 0.08), figure 3(b) illustrates that our approach leads to similar scaling trends. 



For larger values of T+, %nreq,o(^^) deviates from the value obtained using the laminar 



Stokes problem (red dotted curve). We note that for T+ > 150, the solution of Ricco & 
Quadrio ( 2008[ ) provides better agreement with the DNS results. This implies that the 



Reynolds stress that drives the evolution of W in RANS is not as strongly correlated 
with Wy as suggested by the turbulent viscosity of the uncontrolled fiow. In § 4.2 we 



show that accounting for the effect of control on the turbulent viscosity improves the 
prediction of the required power resulting from the use of Wq. 

The apparent lack of infiuence of the wall movements on /7o, observed in § |2.1[ is at odds 
with experiments and simulations that have shown that properly designed oscillations 



can reduce drag by as much as 40% (Jung et a/.||1992 


Baron & 


Quadrio 1996; LaadharT] 


et a/.|1994 Choi et a/.|1998 Choi|2002 Quadrio & Ricco|2003 


2004 


Ricco|2004). Thus, 



model-based control of turbulent fiows requires thorough examination of the infiuence of 
control on vt- For spanwise wall oscillations we address this problem in 



3.2 



2.3. The model equation for vt 

Direct numerical simulations can be used to study the effect of control on turbulent fiows. 
However, resolving all scales of motion at large Reynolds numbers may be prohibitively 
expensive, which motivates use of the Reynolds- averaged equations in conjunction with 
turbulence modeling. The challenge here is to establish a relation between vt and the 
second-order statistics of velocity fiuctuations. By choosing a velocity scale /c^/^ and a 
length scale /c^/^/e, turbulent viscosity can be expressed as ( Pope|[2QQQ ) 

k\y) 



VT{y) 



cRt 



<y) 



(2.12) 



where c = 0.09 is a multiplicative constant (for validity of this assumption in the near- 
wall region we refer the reader to |Pope| ( |2000 )). Both the turbulent kinetic energy k 
and its rate of dissipation e are determined by averaging the second-order statistics of 
fiuctuations over the horizontal directions and one period of oscillations 



Ky) 



1 

2T 
1 
T 



{uu -\- vv -\- ww) (?/,t)dt. 



— / (2 {Uj^Ua; + VyVy + 



ly + WyWy + V^V^ + W , 



+ + v^) {y, t) dt. 



(2.13) 



The most widely used method for computing k and e in engineering fiows is the k-e 
model ( Jones fc Launder]|1972 Launder fc Sharnia||1974 ) where k and e are determined 
by solving two transport equations ( |Pope| 2000 ) . Even though these are less complex 
than the NS equations, they are still computationally expensive, and not convenient for 
control design and optimization. In § |3j we instead develop a simulation- free method. 
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which is computationahy efficient and amenable to control design and optimization, for 
determining the effect of fluctuations on ut in the flow with control. 



3. Stochastically forced flow with control 



Since ut in (2.12) is determined by the second-order statistics of velocity fluctuations, 
we use the stochastically forced linearized NS equations to compute k and e in the flow with 
control. Here, we utilize the fact that the second-order statistics of linear time-periodic 



systems can be obtained from the solution of the corresponding Lyapunov equation ( Jo- 
vanovic fc Fardad||2QQ8 ) . It is well-known that the analysis of the steady- state variance 



of infinitesimal fluctuations around the laminar flow can be used to identify flow struc- 



tures that initiate the onset of turbulence (Farrell & loannou 1993 Bamieh & Dahleh 



2001 Jovanovic & Bamieh 2005). In this paper, we show that eddy- viscosity-enhanced 



linearization of the turbulent flow with control in conjunction with turbulence modeling 
can be used to approximate the turbulent viscosity in a computationally efficient way. 

Next, we examine the effect of control on small- amplitude fluctuations around Uq = 
[ ^o(^) Wo{y^ t)]^ • An equivalent expression for Uq can be found from the steady-state 
solution of the modified NS equations subject to wall-oscillations. 







(U • V) U - VP + {1/Rr) V • ((1 + UTO) (Vu + (Vu)^)) , 



Vu. 



(3.1) 



Model (3.1) is obtained by augmenting the molecular viscosity in the NS equations (2.1) 



with the turbulent viscosity utoj and it facilitates analysis of the dynamics of turbulent 
flow fluctuations (up to first order). 







(Uo • V) V - (V • V) Uo - Vp + {1/Rr) V • ((1 + UTO) (Vv + (Vv)^)) 



(3.2) 

Recent research has demonstrated capability of (3.2) to qualitatively predict the spacing 



and length of near- wall turbulent structures observed in experime nts and simulations (|del 



Alamo & Jimenez 



2006 



by Reynolds & Hussain 



Cossu et al. 



2009' "Pujals et al 2009). This model was used 



(1972) as a method for obtaining closure in the equation for 



velocity fluctuations, and it can be traced back to Townsend (1956). These authors 



noted that the eddy viscosity represents the influence of fluctuations on the background 
turbulent stresses, and that this influence is best captured for the fluctuations that have 
larger wavelengt hs relative to the dominant t urbulent scales. A different interpretation 
was provided by 



del Alamo & Jimenez 



(12006') d 

to model the dissipative effects of small scale turbulent structures on the large scales. 
Our simulation-free design of drag-reducing transverse oscillations involves four steps: 

(i) the turbulent mean velocity in the presence of control is obtained from the RANS 
equations (2.4^ where closure is achieved using the turbulent viscosity of the uncontrolled 
flow 

(ii) k and e are determined from the second- order statistics of fluctuations that are 
obtained from the stochastically forced NS equations linearized around the turbulent mean 
velocity determined in (i); 

(iii) for the flow with control, the modifications to k and e are used to determine the 
modification to the turbulent viscosity, vt] 

(iv) the modified vt is used in the RANS equations 1^2.4] ) to determine the effect of 
fluctuations on the mean velocity, and thereby skin-friction drag and control net effi- 
ciency. 
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Figure 4. Block diagram representing Steps (i)-(iv) of our simulation- free approach for deter- 
mining the effect of control on skin-friction drag in turbulent flows. The hollow arrows indicate 
that some of the coefficients in the model of the output subsystems are determined by the out- 
puts of the corresponding input subsystems. The turbulent mean velocity is updated only once 
in the present study. 



Figure [4] represents these four steps using a block-diagram. We note that the slow 
time evolution of the mean flow (relative to the time evolution of fluctuations) is used to 
separate the update of the mean velocity (Steps (i) and (iv)) from the computation of the 
statistics (Step (ii)) and the update of vt (Step (iii)). Rather than updating the turbulent 
viscosity and the mean velocity in an iterative fashion, we update them only once; in §|4j 
we show that the resulting correction to the mean velocity reliably predicts the optimal 
period of drag-reducing oscillations. Also note that Step (i) amounts to finding the steady- 
state solution of system (3.1) subject to wall oscillations, and that Step (ii) amounts to 



the analysis of the linearized model (3.2) in the presence of stochastic forcing. 



3.1. Computation of the velocity correlations 



The evolution form of the linearized model (3.2 ) that governs the dynamics of fluctuations 
around [UQ{y) WQ{y^t)]^ is given by 



v{y,K,t) 



A{K,t)%l){y,n,t) 
C{n)%l){y,K,t), 



(3.3) 



where i/? = ['^ r^]^ is the state, rj = — is the wall-normal vorticity, and f is the 



stochastic forcing with second-order statistics determined by ( |3.5| ). System ( |3.3[ ) repre- 
sents a A^-parameterized family of PDEs in y and t with time-periodic coefficients. Here, n 
denotes the streamwise and spanwise wavenumbers, k = {kx^ f^z)i and the same symbol 
is used to denote the variables in physical and wavenumber domains (when necessary, we 
will highlight the distinction by explicitly stating the dependence on (x, z) and {kx^ i^z)i 
respectively). The operators A and C in (3.3) are given by 



A = 



All 
A21 





A22 



c = 



1 



IKxdy 














A21 
A22 



{{l/Rr) ((1 
{l/Rr) ((1 + VTO)^ 



(3.4) 



dy) - IK^Uo - IKzWo 
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is the Laplacian, 



and 



where A 

v{±l,n,t) = v'{±l,K,t) = r]{±l,K,t) = 0. 

We next briefly describe a method for determining the steady-state statistics of the 
hnearized system ( |3.3[ ) driven by a zero-mean temporahy white stochastic forcing, with 
second-order statistics, 



(f(-,A^,ti)0f(-,A^,t2)) = M{K)S{ti - t2). 



(3.5) 



Here, S is the Dirac delta function, f f is the tensor product of f with itself, and 
M{n) is a spatial spectral-density of forcing. For homogeneous isotropic turbulence, the 
steady-state velocity correlation tensors can be reproduced by the linearized NS equations 
subject to white-in-time forcing with second-order statistics proportional to the turbulent 



energy spectrum (Jovanovic & Georgiou 2010). Using this analogy, we select M(k) to 



guarantee equivalence between the two-dimensional energy spectra of the uncontrolled 
turbulent flow and the flow governed by the stochastically forced NS equations linearized 
around Uq = [Un( v) 0]^. To this end, we use the DNS-based ener gy spectrum of the 
uncontrolled flow ( del Alamo fc Jimenez||2003 del Alamo et a/. [2004 ), E{y,K,), to define 

EiK) 



M{k) 



Mo{k) 



^E{y,K)I 




Mo{k), 




Here, E{k) = E{y^ k) dy is the two-dimensional energy spectrum of the uncontrolled 
flow, Eo{k,) is the energy spectrum obtained from the linearized NS equations subject to 
a white-in-time forcing f with spatial spectrum Mo(a^), and + denotes the adjoint of an 
operator. Note that the adjoints of the operator appearing in the expression for Mo{k) 
and the operator A should be determined with respect to the inner product that induces 
kinetic energy of flow fluctuations; for additional details, see Jovanovic & Bamieh|( |2005p . 
For the time-periodic system (|3.3[), the operator A in (3.4) can be written as 



A{K,t) = Ao{k) + a{A_i{K)e' 



-iujt t 



Ai{k,) e 



iuJt t\ 



(3.6) 



where the expressions for Aq, and Ai are given in Appendix [B] It is a standard fact 



that the response of the linear time-periodic system (3.3) subject to a stationary input 
is a cyclo-stationary process ( Gardner]! 1990 ), meaning that its statistical properties are 
periodic in time. For example, the auto-correlation operator of ijj is given by 



X{K,,t) = (g)'0(-,K,t)) = 



2 We 



-i2ujt t 



(3.7) 



The averaged effect of forcing (over one period T) is determined by the operator Xq 



hi: 



X{K,t)dt = Xoin). 



(3.8 



In the remainder of the paper, we consider small amplitude of wall oscillations a. 
This choice is motivated by the observation that the power required to maintain the 



oscillations increases quadratically with a (cf. (2.10) and Ricco & Quadrio (2008)). Hence, 
large amplitudes may be prohibitively expensive from control expenditure point of view. 
Furthermore, for sufficiently small value of a the velocity correlations can be computed 
efficiently using perturbation analysis in a ( Jovanovic fc Fardad||2008 Jovanovi6||2008 ) . 
We thus use perturbation analysis to identify wall oscillation periods that yield the largest 
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drag reduction and net efficiency; note that these do not necessarily coincide with each 
other. 



Up to second order in a, the operator Xq in (3.8) is given by 

Xo{n) = Xo,o(k) + a2Xo,2(K) + 0{a^), (3.9) 



where Xq^q and Xo,2 are obtained from the set of decoupled Lyapunov equations (Jo- 
vanovic & Fardad 2008; Jovanovic'2008"); see Appendix [C| for details. The auto-correlation 



operator of the state i/? of stochastically forced uncontrolled flow is determined by Xo,o- 
On the other hand, the operator Xo,2 represents the second-order correction to Xo,o 



induced by wall oscillations. As shown in § 3.2 and § 3.3, Xo,2 determines the effect of 



fluctuations on /c, e, and C/ in the flow with control. 

3.2. Influence of fluctuations on turbulent viscosity 



According to (2.12), ut is determined by the second-order statistics of velocity fluc- 
tuations. By considering dynamics of infinitesimal fluctuations, these statistics can be 
obtained from the auto-correlation operator Xq. In the flow subject to small amplitude 



wall oscillations Xq is given by ( |3.9| ), implying that k and e can be expressed as 

k{y) = ko{y) + aH2{y) + 0{a^), ^ ^ 

(3.10) 

<y) = eo{y) + a^e2{y) + 0{a*). 

Here, the subscript denotes the corresponding quantities in the uncontrolled turbulent 
flow, and the subscript 2 quantifies the influence of fluctuations in the controlled flow at 
the level of a^. A computationally efficient method for determining k2 and 62 from Xo,2 
is provided in Appendix [P) 

For small amplitude oscillations, substituting k and e from ( 3.10[ ) into ( |2.12 ) yields 



. . 02 fc!M = p2 (fco(y) + a'hiy) + Oia^)f 
vtKV) cti, ^^^^ CH, ^^^^^ ^ ^ , 

which in conjunction with Neumann series expansion leads to 

^2k2{y) e2{y)\ (3.11) 



\ko{y) eo{y)^ 

Therefore, up to second order in a, the influence of fluctuations on turbulent viscosity in 
the flow with control is determined by second-order corrections to the kinetic energy /c2 
and its rate of dissipation 62- 

3.3. Skin-friction drag coefficient and net efficiency 

We next show how velocity fluctuations in the flow subject to small amplitude oscillations 
modify the skin-friction drag coefficient and the net efficiency. As discussed in § |2.2[ C/ is 
determined by U and %nnet is determined by both U and W. The influence of fluctuations 
on U and W in the flow with control can be obtained by substituting ut from ( |3.11[ ) 



into (2.4), and thereby expressing /7, W, and Px as 

U{y) = Uo{y) + a^U2{y) + 0{a^), 
Wp{y) = WM + o^^WpAy) + 0{a^)^ (3.12) 
P. = -1 + + 0{a^), 

where the expressions for ^p,2^ and Px^2 are provided in Appendix IE] 
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An expression for the saved power is obtained by substituting Px from (3.12) into (2.9) 

%Cf = %nsave = C.2%nsave,2 + %nsave,2 = 100P,,2. 

In flows subject to small amplitude oscillations, the above equation shows that a positive 
(negative) value of %nsave,2 signifles drag reduction (increase). On the other hand, the 
required power can be obtained by substituting Wp from (|3.12) into (2.10) 



%nreq = 



req,0 



%n 



req,2 ) 



(3.13) 



where %nreq,o is given in (2.11), and %nreq,2 is provided in Appendix [e} We note the 
contrast in the way fluctuations influence saved and required powers; while fluctuations 
make O(a^) contribution to %nsave5 they affect %nreq only at the level of a^. This 



explains the close agreement, observed in 



2.2 



between %nreq,o and the DNS results 



of Quadrio & Ricco (2004). Finally, the net efficiency is given by 

%nnet = C.'%nnet,2 + 0{a^), %U^et,2 = %nsave,2 - %nreq,0. 

The developments of this section are used in § |4] to determine the skin-friction drag 
coefficient and the net efficiency in the ffow subject to wall oscillations. 



4. Turbulent drag reduction 

In this section, we examine the effect of transverse wall oscillations on skin- friction 
drag and net efficiency in ffows with R^- = 186, 547, and 934. For these Reynolds num- 
bers, the second-order statistics of the uncontrolled turb ulent ffow were obt aine d using 



3.1 



we use 



DNS ddel Alamo fc Jimenez| |2003 del Alamo et al. ] |2004| ). As explained in 
available DNS data t o de termine the sp atial spectrum of the stochast i c forcing (|3.5[) in 
the evolution model (3.3). The results of 



(2003); 



del Alamo et al. 



' del Alamo & Jimenez 
(|2004|) also provide the turbulent kinetic energy in the uncontrolled flow, /cq, and thereby 
its rate of dissipation. 



eo(^) 



cR 



where uto is deffned by (2.8) 



The differential operators in the wall-normal direction are discretized using Ny col- 
location points (Weideman & Reddy 2000). In horizontal directions, we use Nx x A^^ 

, where k,x 



wave numbers with < a^^^ < tvx,ma^ and < 
the largest wavenumbers used in the DNS of 



and 



del Alamo & Jimenez 



(2003); 



del Alamo 



et al. (2004); table [T] provides summary of parameters used in our study. The DNS-based 
energy spectrum taken from http://torroja.dmt.upm.es/ftp/channels/data/ is interpolated 
on the A^^ collocation points in y. 

4.1. Saved power 

We ffrst examine the effect of period of oscillations T+ on the turbulent drag reduction 
and the saved power. Solid curve in figure [5] shows the second-order correction to the 
saved power %nsave,2(^^)7 normalized by its maximum value max (%nsave,2(^^)) = 2.8, 
for the controlled fiow with R^ = 186. The positive value of %nsave,2 indicates that drag 
is reduced for all values of that we considered, with the largest drag reduction taking 
place at T+ = 102.5. Our theore tical predictions are in close agreement with DNS at 
Rr = 200 ( [Quadrio fc Ricco||2004 ), where it was shown that T+ ^ 100 yields the largest 
drag reduction for control amplitudes a = 3.1 and 6.2. 

Figure [5] also compares %nsave,2 (solid curve) obtained using our analysis with %nsave 
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Rr 



Cl 



C2 



Ub 



186 
547 
934 



50 
50 
50 



101 
151 
201 



51 
51 
51 



42.5 
128 
255 



84.5 
255.5 
511.3 



46.2 
29.4 
27.0 



0.61 
0.45 
0.43 



15.73 
18.38 
19.86 



Table 1 . The parameters used in our study. At each Rr , ci and C2 are select ed to mini mize the 
least squares deviation between t he mean streamwise velocity obtained fr om (|2.5|)-(|2.8|) and the 
mean velocity obtained in DNS (*del Alamo Jimenez"2003| 



flux t/s is kept constant by adjusting the pressure gradient. 



del Alamo et a/.||2004p . The bulk 




Figure 5. (Solid curve) Second-order correction to the saved power, %nsave,2(T^) normalized 
by max (%nsave,2(T^)) = 2.8, as a function of the period of oscillations for the flow with 
Rt = 186; (Symbols) DNS-based % nsave(T^) normalized by the corresponding largest values at 
= 200 ( iQuadrio Ricco|20"04| for control amplitudes a = 2.25, o; a = 6, □; and a = 9, V. 



(symbols) obtained using DNS at Rr = 200 ( Quadrio fc Ricco|2004 ); both sets of results 
are normalized by their corresponding maximal values. In DNS, the largest drag reduction 
takes place at T+ ^ 100 for a = 2.25 and a = 6, and at T+ ^ 125 for a = 9. Thus, 
for small control amplitudes, perturbation analysis up to second order in a reliably 
predicts optimal period of drag reducing oscillations. For the optimal period of oscillations 
and a = 2.25, our perturbation analysis predicts 13.6% drag reduction, whereas 17.4% 
drag reduction is obtained in DNS (Quad rio fc Ricco||2QQ4 ). The quantitative difference 
between the DNS results and the results of perturbation analysis may be attributed to 
the effects of higher order corrections. Another factor that warrants further scrutiny is 
modeling of the spatial spectrum of stochastic forcing. Analysis of these effects is beyond 
the scope of the current study. 

For small oscillation periods, T+ ^0.5, the second-order correction to the saved power 
is negligible. Recent research has determined necessar y conditions on the amplitude and 
period of oscillations for drag reduction. For T+ > 30, Ricco fc Quadrio] pQQS ) used the 



solution to the Stokes second problem to show that the smallest oscillation amplitude that 
yields drag reduction is approximately equal to 1.8 (in viscous units). It was also shown 
that the amplitudes of oscillations that achieve drag reduction become prohibitively large 
as the Stokes layer thickness 5s ~ a/T+ approaches zero. In addition, Quadrio fc Ricco 
(2011) showed that the smallest value of the Stokes layer thickness that results in drag 
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Figure 6. (Solid curve) Fourth-order correction to the required power, %nreq,2(T^), as a 
function of the period of oscihations for the flow with Rt = 186; (Filled symbol s) the 
difference between the DNS-based required power %nreq (shown by o pen symbols in figure js^aj) 
and the required power %nreq,o (shown by solid curve in figure |3 (a) |) that is obtained from the 
solution Wo (y, t), %U 

req 

Ricco|f2004| for a = 

max(|7ollreq,2(- 



- a"%nreq,o(T+). %nreq Is obtalucd usmgT)NS at R 
2.25, •; a = 6, ■; and a = 9, T. 

req,2tT+)|) = 0.0126 and the largest values of |%nreq 



200 dQuadrio 
The results are normalized 
-a2%n,eq,o(T+)|. 



reduction by traveling waves of spanwise wall motion is approximately equal to 1 in 
viscous units. 

4.2. Required control power 



We next study the power required to maintain wall oscillations. From (3.13) it follows 



that, up to second order in a, %nreq,o defined by (2.11) determines the required power, 
we also recall that the discrepancy between %nreq,o 



From 



2.2 



and the DNS-based 



required power gets larger for > 150; cf. figure 3(a) Here, we show that accounting 



for the effect of fluctuations in the flow with control reduces this discrepancy. As shown 
in figure |6j the fourth-order correction to the required power, %nreq,2(^^), is negative 
for > 45 and it decreases with T+. This is in agreement with our earlier observation 
that %nreq,o Overestimates the required power obtained in DNS. The filled symbols in 



figure [6] represent the difference between the required power obtained by |Quadrio fc Ricco 



%nreq 

power, %n 



(2004) using DNS (shown by open symbols in figure 3(a)) and the required power that 
we obtain using %nreq,o (shown by solid curve in figure 3(a)). For a = 2.25, 6, and 9, 
a^%nreq,o IS in closc agreement with the fourth-order correction to the required 
req, 2 (soHd curvc in figure |6|. 

4.3. Net efficiency 

The net efficiency in the flow subject to wall oscillations is determined by the difference 
between the saved and required powers. The solid curve in figure [7| shows the second- 
order correction to the net efficiency, %nnet,2(^^), in the flow with = 186. We see 
that %nnet,2 > for T+ > 75, indicating that, for small control amplitudes, a positive 
net efficiency can be achieved if the period of oscillations is large enough. Our prediction 
is in close agreement with DNS at Rr = 200 (Quadrio fc Ricco 2004|) where positive 
net efficiency of oscillations with a = 2.25 is obtained for T+ > 70. Furthermore, up to 
second order in a, the largest net efficiency takes place at T+ = 180. This value differs 



from the value of T+ that yields the largest saved power, T+ = 102.5 (cf. § 4.1). This 
difference can be explained by the fact that the peak of %nnet,2 = %nsave,2 — %^req,o^ 
takes place at a value of T+ where the derivatives (with respect to T+) of %nsave,2 and 
%nreq,o(T+) are equal to each other. Since %nreq,o(^^) is a monotonically decreasing 
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%nnet,2(T+);%nnet/a' 




280 



Figure 7. Comparison between second-order correction to the net efficiency %nnet,2(T^) (soiid 
curve) for ttie ffow witfi Rr — 186 (soiid c urve), and %Ilnet/a'^ (symbois). Symbois stiow DNS 
data at Rr = 200 ( |Quadrio Ricco||2004t for a = 2.25, o. 



function, its derivative is negative for all T+. Therefore, the largest net efficiency takes 
place at some T+ > 102.5 where the slope of the curve %nsave,2(^^) is also negative. 
(Our analysis shows that this happens at T+ = 180.) The symbo ls in figure [t] show 
%nnet/Q^^ obtained from DNS at Rr = 200 ( [Quadrio fc Ricco||2004[ ) for a = 2.25. Even 
though the essential trends are captured by %nnet,2 (solid curve), the DNS net efficiency 
peaks at T+ = 125, which is approximately 30% smaller than the value of T+ predicted 
by our perturbation analysis. This discrepancy may be attributed to a slower rate of 
decay of %nsave,2 relative to %nsave obtained in DNS for T+ > 100; cf. figure [sj 

We note that in DNS the net efficiency becomes negative for large control amplitudes. 
At Rr = 200, Quadrio & Ricco (2004) showed that %nnet becomes negative for all T+ if 



a > 3.5. Therefore, the positive net efficiency predicted by perturbation analysis is only 
valid for small control amplitudes. This can be explained by noting that perturbation 
analysis predicts quadratic increase of both saved and required powers with a. On the 



other hand, Quadrio & Ricco (2004) showed that saved power exhibits slower than linear 



growth with a for large control amplitudes. 

4.4. Drag reduction in flows with larger Reynolds numbers 

After assessing utility of perturbation analysis for flows with Rr = 186, we turn our 
attention to the effect of control at Rr = 547 and 934. Figure |8(a)| shows that the 
second-order correction to the saved power, %nsave,2, is positive for all Rr, and that the 
optimal T+ slightly decreases with Rr (T+ = 102.5, 96, and 94 for Rr = 186, 547, and 
934, respectively). The corresponding periods of oscillations in outer units, T = /Rr, 
are T = 0.55, 0.18, and 0.10. Therefore, as Rr increases, larger frequency of oscillations is 
required for the optimal drag reduction. In addition, since %nsave,2 decreases with Rr, the 
drag-reducing ability of wall oscillations deteriorates at higher Reynolds numbers, which 
is in agreement with the observations of Choi et al. (2002); Ricco & Wu (2004); Ricco 



fc Quadrio (2008); Touber fc Leschziner (2012). Our results suggest that the amount 
of drag reduction scales as R~^-^^, which a grees fairly well with the d ecline in drag 
reduction R~^ '^ found by Choi et al. (2002); Touber fc Leschziner (2012). Our analysis 
also demonstrates that for T+ = {40, 100, 160, 220} the second order correction to the 
saved power reduces by {10.4%, 20.2%, 20.8%, 18.8%} when Rr increases from 186 to 
547, and by {16.6%, 27.2%, 28.1%, and 26.1%} when Rr increases from 186 to 934, 
respectively. This is in line with deterioration of 18% observed in the DNS of Touber 
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Figure 8. The second-order correction to (a) the saved power, %nsave,2(T^); (b) the required 
power, %nreq,o(7'^); and (d) the net efficiency, %nnet,2(T^); and (c) the fourth-order correction 
to the required power, %nreq,2(T^), as a function of the period of oscihations for the flows 
with Rr = 186, 547, and 934.' 



& Leschziner (2012) when increases from 200 to 500. In addition, our predictions 
qualitatively agree with the results of Ricco & Quadrio ( 2008[ ) where it was shown that 
deteriorating effect of Rr on drag reduction becomes stronger with increase in T+. 

Figure 8(b) shows monotonic decrease of %nreq,o(^^) with both T+ and R^-. We also 
note that a product between the bulk flux Ub and %nreq,o does not change with R^. This 
demonstrates invariance under change in R^ of the second-order correction to the required 
power (before normalization is done). [Ricco fc Quadrio" ( 2008| ) noted that since Ub 
approximately scales as R^-^^^^ the required power should approximately scale as R~^-^^^. 
For T+ = {40, 100, 160, 220} and r = {0.145, 0.141, 0.138, 0.135}, respectively, we 
indeed obtain very weak variations of %nreq,o/^r ^ with R^. The fourth-order correction 
to the required power, %nreq,2 5 is also reduced as Rr increases; see figure 8(c)| The 
difference between the quantities shown in figures |8(a)| and 
order correction to the net efficiency, %nnet,2- From figure 
net efficiency reduces with Rr, and that %nnet,2 saturates for large Reynolds numbers 



8(b) 



8(d) 



determines the second- 
we see that the largest 



20 R. Moarref & M. R. Jovanovic 




Figure 9. (a) The turbulent viscosity, z^to(2/^); (b) the turbulent mean streamwise 
velocity, Uo{y'^), in the uncontrolled flows with Rr = 186, 547, and 934. 

i^T2{y+;T+) U2{y+;T+) 




Figure 10. Second-order correction to (a) the turbulent viscosity, VT2{y^ ]T^)] and (b) the 
mean streamwise velocity, U2{y~^;T~^), for Rr = 186 and different values of 30 < < 300, 
where increases in the direction of the arrows. The thick curves correspond to the that 
yields the largest drag reduction (cf. figure 11). 



4.5. Effect of control on turbulent viscosity and turbulent mean velocity 

We next examine the effect of wall oscillations on the turbulent viscosity and the turbulent 
mean velocity. Figure [9] shows lyroiv^) and Uo{y^) for Rr = 186, 547, and 934. We note 
that the profiles ofUo for different Rr lie on the top of each other, and that uto does not 
scale in wall units. In particular, the peak of uto takes place at y~^ ^ Rr/"^^ approximately 
half way between the walls and the channel centerline. On the other hand, the effect of 



control on ut is strongest in the viscous wall region, y~^ < 50; see figure 10(a) 



Figure [10] shows the second-order corrections to ut2 and U2 for Rr = 186 and 30 < 
T+ < 300. Since ut2 < for all T+, perturbation analysis up to second order in a 
predicts turbulence suppression for all periods of oscillations; see figure |10(a)[ Further- 
more, turbulence suppression region shifts away from the walls with increase in T+. 
We note that su ppression of the turbulen t burst ing process was observed in DNS for 
25 < T+ < 200 dJung et a/.]|l992[ ). Figure 1 10 (b)] shows that the oscillations reduce the 
mean velocity gradient in the immediate vicinity of the walls {U2 < for y^ < 13). 



Model-based design of transverse wall oscillations for turbulent drag reduction 21 

^T2(?/+) U2{y^) 




Figure 11. The second-order correction to (a) the turbulent viscosity, h'T2{y^)] and (b) the 
mean streamwise velocity, U2{y'^), for the values of that yield the largest drag reduction in 
the flows with Rr = 186, T+ = 102.5; Rr = 547, T+ = 96; and Rr = 934, T+ = 94. 



On the other hand, the mean velocity is shifted upward in the log-law region. Both these 



trends were previously observed in the experiments (Choi 2002) and DNS (Choi et al. 
1998 Baron fc Quadrio]|1996 ). Choi (2002) argued that the negative spanwise vorticity 



introduces near wall modifications to the mean streamwise velocity, and that it suppresses 
the production of turbulence by weakening the vortex stretching mechanism. Our per- 
turbation analysis underestimates t he value of above which the shift-up in U takes 
place; y+ > 13 vs y+ > 30 in DNS jBaron fc Quadrio||l996| [Tbuber fc Leschziner||2012 ) 
and experiments ( |Choi et a/.||1998[ [Laadhari et a/.||1994[ |Ricco fc Wu||2004D . 



Figure 



11 



compares U2 and z/t2 at three Reynolds numbers for the values of T+ that 
induce the largest drag reduction. We see that z^t2 peaks at y~^ ^ 20 for all Rr. This 
suggests that the optimal drag-reducing frequency minimizes the turbulent viscosity near 
the interface of the buffer layer and the log-law region. Even though the negative peak of 
ut2 increases with R^- (cf. figure 11(a) ), the ratio of ut2 and uto decreases with R^. Thus, 
wall oscillations are less effective in suppressing turbulence at larger Reynolds numbers. 
Finally, figure pT(b)| shows that the slope of U2 decreases with R^-] this is in agreement 
with our earlier observation that smaller drag reduction is achieved at higher Reynolds 
numbers. 



4.6. Effect of control on turbulent kinetic energy 



We next examine the effect of control on the kinetic energy of fluctuations. Figure 12 



compares the premultiplied two-dimensional energy spectrum of the uncontrolled flow, 
nxi^zE{K,)^ with the premultiplied second-order correction to the energy spectrum, K,xf^z 
E2{i^), in the flow subject to wall oscillations with the optimal drag-reducing period 
T+ = 102.5. The energy spectra are premultiplied by the spatial wavenumbers such 
that the area under the log-log plot is equal to the total energy of fluctuations. We 
see that the most energetic modes of the uncontrolled flow take place at hZx ^ 2.5 and 
Hiz ^ 6.5; cf. figure 12(a)[ In addition, wall oscillations further amplify the most energetic 
modes of the uncontrolled flow with small streamwise wavelengths (cf. red regions in 
figure [12(b)]) , while they suppress the most energetic modes of the uncontrolled flow with 
large streamwise wavelengths (cf. blue regions in figure [T2(b)] ). This agrees well with the 
study of |Ricco| ( |2011 ) that examined the effect of spanwise wall oscillations on free-stream 
vortical structures in a Blasius boundary layer. Over a range of spanwise wavelengths. 
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Figure 12. (Color online) ( a) Premultiplied DNS-base d energy spectrum of the uncontrolled 
flow, k,xI^zE{k), at Rr = 186 (del Alamo & Jimenez 2003 ); and (b) second-order correction to the 
energy spectrum, hixi^zE2{K,), for the How subject to wall oscillations with optimal drag-reducing 
period T+ = 102.5. 

it was found that wall oscillations amplify (attenuate) disturbances with small (large) 
streamwise wavelengths. Figure [T2(b)| shows that the largest energy amplification takes 
place at hCx ^ 4.4 and ^ 10.2, and the largest energy suppression occurs at hix ^ 0.8 
and ^ 8- 8- The total effect of control on the kinetic energy can be quantified by 
Ji^E2{k,) dn/ Ji^E{k,) dn^ which for wall oscillations with T+ = 102.5 is approximately 
— 1.5%. This yields 7.5% reduction in the total energy of fluctuations for a = 2.25. 

The effect of wall oscillations on the turbulent kinetic energy and its rate of dissipation 
is shown in figure 13 We see that the second-order correction to the kinetic energy, /c2, 
is negative for T+ = 102.5, suggesting that the turbulent kinetic energy is reduced 
with the largest suppression taking place at ~ 8.6; cf. figure 13(a) In addition. 



the second-order correction to the rate of dissipation of turbulent kinetic energy, 62, is 
negative almost everywhere (except for a small region 5.4 ^ ^ 8.9) and the largest 
reduction occurs in the viscous sublayer < 5; cf. figure 13(c) Therefore, perturbation 
analysis up to second order in a captures previously made experimental and numerical 
observations that wall oscillations suppress both the production and dissipation of the 
turbulent kinetic energy ( [Jung et a/.|1992 Choi et a/.|1998 Baron fc Quadrio|1996 Choi 
2002 ). Figures |13(b)| and |13(d~ show that, relative to the uncontrolled flow, the wall 
oscillations with T+ = 102.5 and a = 2.25 have more profound influence on the turbulent 
kinetic energy than on its rate of dissipation. This suggests that the turbulent production 
is suppressed more than the turbulent dissipation, which explains the reduced turbulent 



viscosity ut in the flow with control; cf. figure 11(a) We observe close agreement between 
the modification to the turbulent kinetic energy in figure |13(b)| and the modification to 
the dominant component of the Reynolds normal stress reported in the DNS of Touber 
fc Leschzinerl (|2QT2|) at Rr = 500 (their figure 5(b)). 



5. Turbulent flow structures 



In this section, we use stochastically forced linearized model (3.2) to examine the effect 
of wall oscillations on the turbulent flow structures. We only present results for i^^- = 186 
and note that similar flow structures are observed for all Reynolds numbers that we have 
considered. 
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Figure 13. (Color online) (Left column) Second-order correction to (a) the turbulent kinetic 
energy k2{y~^), and (c) its rate of dissipation 62(^^)5 for the flow subject to wall oscillations 
with optimal drag-reducing period = 102.5 at Rr = 186. (Right column) (b) Comparison 
between the turbulent kinetic energy in the uncontrolled flow ko (black) and in the flow with 
control ko + a^k2 (dashed red); (d) Comparison between the rate of dissipation of turbulent 
kinetic energy in the uncontrolled flow eo (black), and in the flow with control eo + 0.^62 (dashed 
red), for a = 2.25. 



We decompose the velocity field into characteristic eddies (Moin & Moser 1989) by 
determining the spatial structure of the fluctuations that contribute most to the kinetic 
energy at a given k = {f^x^f^z)'', see Appendix [F| for details. It is worth noting that the 
dominant characteristic eddy resulting from the analysis of the stochastically forced lin- 
earized model in the uncontrolled flow qualitatively agrees with the results obtained using 



eigenvalue decomposition of the DNS-based autocorrelation matrices; compare figure [15 
of this paper with figure 15 in Moin & Moser (1989). 

We next examine the effect of wall oscillations on the dominant characteristic eddy. 
Figure 14 compares the streamwise streaks and their surrounding vortex core for the 
dominant characteristic eddies in the uncontrolled flow (left column) and in the flow 
subject to wall oscillations with T+ = 102.5 and a = 2.25 (right column). In both cases, 
the iso-surfaces represent high- (red) and low- (blue) speed streaks at 70% and 60% of 
their largest absolute values in the uncontrolled flow, respectively. The vortex core (green 
surface) is obtained based on the 'swirling strength' criterion which identifies motions 



with large rate of rotation and large orbital compactness ( Chakraborty, Balachandar 
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Figure 14. (Color online) Three-dimensional iso-surfaces of the streamwise streaks (red and 
blue) and the vortex core (green) for the characteristic eddy in the uncontrolled flow (left 
column) and the flow subject to wall oscillations (right column) with optimal drag-reducing 
period = 102.5 at Rt = 186. (a)-(b): bird's-eye view; (c)-(d): side view; (e)-(f): front view; 
(g)-(h): top view. The fast- (red) and slow- (blue) moving streaks are respectively shown at 70% 
and 60% of their largest values in the uncontrolled flow, and the vortex core is obtained based 
on the 'swirling strength' criterion ( |Chakraborty et a/.|2005| ) with Xd > 12, \Xcr/Xci\ < 0.4, and 

Acz/12- |Acr/Ac^|/0.4= 1.2. 
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Figure 15. (Color online) Cross sections of the streamwise velocity (colored contours) and 
the spanwise vorticity (black contours) for the characteristic eddy in the uncontrolled flow (left 
column) and the flow subject to wall oscillations (right column) with optimal drag-reducing 
period T+ = 102.5 at Rr = 186. (a)-(b): z+ = 0; (c)-(d): x+ = 0; (e)-(f): = 3.8; (g)-(h): 
= 10.8. The streamwise velocity is normalized by its largest absolute value in the uncontrolled 
flow, and the level sets for the spanwise vorticity correspond to 40%, 60%, and 80% of its largest 
value in the uncontrolled flow. 
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Adrian|[2QQ5 ) . This criterion can be expressed in terms of the real and imaginary parts 



of the complex eigenvalues, Xcr iiAci, of the rate of strain tensor at each point inside the 
channel. The large rate of rotation requires Xd > 61, where bi determines the strength 
of the swirling motion. In addition, orbital compactness is guaranteed if \Xcr/Xci\ < ^2, 
where 62 determines the compactness of the swirling motion in the plane spanned by the 
real and imaginary parts of the eigenvector corresponding to Xcr i iAci- For example, 
62 = identifies a pure circular motion, and larger values of 62 allow for inclusion of the 
in-plane converging or diverging spiral motions to the vortex core. In figure 14, we use 



bi = 12 and 62 = 0.4 to identify the strong vortex core that surrounds the slow-moving 
streaks. We see that wall oscillations reduce intensity and spatial spread of both the 
streamwise streaks and the vortex core, which is in agreement with experiments of .Ricco] 
( |2QQ4| . 

Figure 15 shows the streamwise velocity, u{x^y^z)^ and the spanwise vorticity, uj^ = 
Vx{x^ z) — Uy{x, z), for the characteristic eddy in the uncontrolled fiow (left column) 
and in the flow subject to wall oscillations with T+ = 102.5 and a = 2.25 (right column). 
The spanwise vortices are shown at 40%, 60%, and 80% of their largest values (black 
contours). The streamwise streaks (colored contours) are normalized by their largest ab- 
solute value. We see that, in the viscous sublayer < 5, wall oscillations suppress the 
largest spanwise vorticity by approximately 6.5%. In addition, the magnitude of stream- 
wise streaks is reduced by approximately 12%. Choi (2002) argued that transverse wall 
movements induce negative spanwise vorticity in the flow, thereby suppressing intensity 
of the streamwise streaks. In addition, wall oscillations slightly move the center of the 
spanwise vortices away from the wall {y~^ = 3.8 vs. y~^ = 4.3), and they signiflcantly 
reduce the upstream extent of the spanwise vortices. 



6. Concluding remarks 

This paper has introduced a model-based approach to controlling turbulent flows. In 
contrast to standard practice that embeds turbulence models in numerical simulations, we 
have developed a simulation-free approach that enables computationally-eflicient control 
design and optimization. This has been achieved by merging turbulence modeling with 
techniques from linear systems theory. In particular, we have used the turbulent viscosity 
hypothesis in conjunction with a model equation for ut based on k and e to determine 
the influence of turbulent fluctuations on the mean velocity in the flow with control. 

We have shown that the study of dynamics is of prime importance in designing drag- 
reducing wall oscillations. This has allowed us to determine the influence of control on 
the turbulent viscosity in a simulation-free manner. This contribution goes well beyond 
the problem that was used to demonstrate the predictive power of our model-based 
control design - turbulent drag reduction by transverse wall oscillations. The computa- 
tional complexity of determining the turbulent viscosity (in the flow with control) has 
been signiflcantly reduced by obtaining k and e from the second-order statistics of eddy- 
viscosity-enhanced linearized model with stochastic forcing. 

The flrst step in our control-oriented modeling involves augmentation of the molecular 
viscosity with the turbulent viscosity of the uncontrolled flow. The resulting model is then 
used to determine the turbulent mean velocity in the flow with control, and to study the 
dynamics of velocity fluctuations around it. By considering linearized equations in the 
presence of white-in-time stochastic forcing (whose spatial spectrum is selected to be 
proportional to the turbulent kinetic energy of the uncontrolled flow), we have quantifled 
the influence of control on the second-order statistics of velocity fluctuations and thereby 
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on the turbulent viscosity. Finally, the modifications to the turbulent viscosity determine 
the turbulent mean velocity and skin-friction drag in the flow with control. 

Since the evolution model for flows subject to wall oscillations is time periodic, a wide- 
sense stationary stochastic forcing induces velocity fluctuations with cyclo-stationary 
statistics. Computing these statistics is challenging even in the linearized case. Motivated 
by the observation that large control amplitudes yield poor net efliciency, we have used 
perturbation analysis (in the amplitude of oscillations) to quantify the effect of control 
on the turbulent statistics in a computationally eflicient manner. 

We have shown that perturbation analysis up to second order reliably predicts the op- 
timal period of drag-reducing oscillations. Furthermore, even though the required power 
obtained using the turbulent viscosity of the uncontrolled flow agrees well with the val- 
ues obtained in DNS, this agreement has been further improved by accounting for the 
effect of control on fluctuations (and, consequently, on turbulent viscosity). In addition, 
the predicted net efficiency resulting from perturbation analysis qualitatively agrees with 
the DNS results. Perturbation analysis has also captured suppression of the turbulent 
kinetic energy and its rate of dissipation by wall oscillations, as well as modiflcations 
to the streamwise component of the turbulent mean velocity (reduction in the viscous 
sublayer and buffer layer and increase in the log-law region). Finally, the spatial spectral 
density tensors of velocity fluctuations obtained from the solution of the corresponding 
Lyapunov equations determine the effect of control on the dominant flow structures. As 
previously observed in experiments, and conflrmed by our analysis, wall oscillations re- 
duce the spatial spread and magnitude of the dominant characteristic eddies and suppress 
the spanwise vorticity in the viscous sublayer. 

It is noteworthy that simple turbulence modeling (that relies on the turbulent viscosity 
hypothesis with expressed in terms of k and e) in conjunction with eddy- viscosity- 
enhanced linearization of the flow with control has significant predictive power for captur- 
ing full-scale phenomena. Even though this model does not reveal all aspects of turbulent 
fiow physics, we have shown that it is well-suited for control design and optimization. 
Development of more sophisticated control-oriented turbulence models may further re- 
duce the gap between theoretical predictions and experiments/simulations. In addition, 
the predictive power of the proposed approach can be enhanced by optimization of the 
power spectrum of the forcing. We expect that our model-based approach will find use 
in designing feedback-based and sensor-less turbulence suppression strategies in other 
geometries, including pipes and boundary layers. 



Acknowledgments 

Financial support from the National Science Foundation under CAREER Award CMMI- 
06-44793 and from the University of Minnesota Initiative for Renewable Energy and the 
Environment under Early Career Award RC-0014-11 is gratefully acknowledged. The 
University of Minnesota Supercomputing Institute is acknowledged for providing com- 
puting resources. The authors would also like to thank anonymous reviewers and the 
associate editor for their valuable comments. 



Appendix A. The required power to maintain wall oscillations 

The wall oscillations require an input power to balance the spanwise shear stresses at 
the walls. The required power over one period T per unit area of the channel walls is 
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obtained from ( Currie||2QQ3 ) 



req 



1 

T 



{W{y,t)r23{y,t))dt 



2/ = ! 



1 

T 



{W{y,t)r23{y,t))dt 



where r23(^,t) = /j^W^y^t) denotes the spanwise shear stress and fi is visc osity . An 
equation for Ilreq (normahzed by pu'^) can be obtained by substituting W from (2.6) and 
using the boundary conditions on Wp given in (2.7) 



n -2"' 

Alreq — 7-, 



P\y=l 



Relative to the power necessary for driving the uncontrohed flow, the required power is 
given by %nreq = 100 nreq/(2/7B), which yields (|2.10|). 



Appendix B. Operators Aq, Ai, and A_i in (3.6) 



The operators Aq, Ai, and A-i in (3.6) are obtained by substituting W from (2.6) 



into the expression ( |3.4| ) for A, 

Ao = 



^0,11 



, ^0,22 = (1/i^r) ((1 + i^To)A + u^^dy) - m^^Uo, 



i/.,A-i (iy;;;^o - ^m,oA) o 



Here, VFm,o(^) = ^^,0(^)5 ^o, ^^,0, ^to and their derivatives denote multiplication 
operators in the wall-normal direction. 



Appendix C. Computing the velocity correlations 

For the time-periodic system (3.3 ), the normal modes are determined by Bloch waves (|^eh 
&: Keller]|1964 Bensoussan, Lions h Papanicolaou||1978 ) 

i{y,K,t) = ^f„(2/,«)e'(^+"-*)*, 



n G Z 



where G [0,cj^) is the angular frequency. The frequency response of the time-periodic 
system ( |3.3[ ) is an operator that maps the bi-infinite input column vector col {fnjnez, into 
the bi-infinite output column vector col{v^}^^^ (Jovanovic & Fardad 2008[ |Jovanovi6 



2008). The system states can also be defined as a bi-infinite column vector colji/j^jnez- 
As discussed in § 3.1 • , a^, t) is a cyclo-stationary process with second-order statistics 
given by (3.7). The kernel representation Kx^{y,£,, t^) of the auto-correlation operator 
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Xr of • , A^, t), 



{il){y,K,t)il)''{^,K,t)) 



E 



i rojt t 



can be expressed in terms of {ijj^}nez 



E 

n G Z 



Furthermore, the frequency representation of the auto-correlation operator of ,K,^t) 
is a self- adjoint bi-infinite block- Toeplitz operator that is parameterized by k 



XiK) = Toep{...,X+,X+,|Xo|,Xi,X2,...}, 

where the box denotes the element on the main block diagonal of A'. 

For the case where f is a zero-mean white process in y and t with second-order statistics 
given by (3.5), we have fo = f, and = for n 7^ 0. Thus, the frequency representa- 
tion of the spectrum of f is given by a bi-infinite block-diagonal operator A4{k) with 
block diagonals equal to M{k). The auto-correlation operator of the state, A'(k), can 
be obtained from the following Lyapunov equation ( Jovanovic fc Fardad|2QQ8[ | Jovanovi6 



2008) 



Here, ^ is a bi-infinite block-diagonal operator 

g{0) = diag{i (i9 + nujt)l}nez. 
and ^ is a bi-infinite block- Toeplitz operator 

A{k,) = Toep {. . . , 0, aAi,\ Aq |, 0, . . .}. 



(CI) 



The solution to (|C l|) can be efficiently computed using perturbation analysis in a (Jo- 



vanovic & Fardad 2008; Jovanovic 2008). The operator T is decomposed into a block- 



diagonal operator J^q and an operator J^i that contains the first upper and lower block 
sub-diagonals 



^0 
^1 



diag{Ao - inutl} 



n G ^ 



(C2) 



Toep{...,0,Ai, ,A_i,0,...}. 



For sufficiently small a, the solution to (CI) can be written as (Jovanovic & Fardad 
2008] ) 

A' = A'o + aA'i + A's + + .... (C3) 

Substituting ( C 2 ) and ( C 3 ) into (CI) and collecting equal powers of a yields the fol- 
lowing set of Lyapunov equations 

Since is block-diagonal, A'^ inherits the structure of the right-hand-side of the equa- 
tion at 0{a^). The structure of the above equations reveals that is a self-adjoint 
block-diagonal operator, is a self- adjoint block- Toeplitz operator where only the first 
upper and lower block sub-diagonals are non-zero, and A2 is a self- adjoint block- Toeplitz 
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operator where only the main block diagonal and the second upper and lower block 
sub-diagonals are non-zero 



^eo{f^) = Toep{ 
A'^i(a^) = Toep{ 
X02{n) = Toep{..., 







.., , 

, 



-^0,0 



^2,2' 



X, 



0,2 



, 
, 



, 

-^2,2, 







...}. 



The above structure of the operator A' in conjunction with the fact that only the element 
on the main block diagonal of A' contributes to the averaged effect of forcing on the 
velocity correlations (cf. equation (3.8)) reveal that, up to second order in a, only Xo,o 
and Xo,2 contribute to Xq 

Xo{k) = Xo,o(a^) + a'Xo,2(A^) + 0{a^), 

The operators Xo,o and Xo,2 are obtained from a set of decoupled Lyapunov equations 



whose size is equal to the size of each block in the bi-infinite Lyapunov equation ( C 1 ) ( Jo- 
vanovic fc Fardad||2008D 



Ao Xo,o + Xo,o = -M, 
(Ao + ic^t/)Xi,i + Xi,iA+ = -(A_iXo,o + Xo,oA+) 

^0 



AoXo,2 + Xo,2A+ = -(A_iX+i + Xi,iAti + AiXi,i + X+i 



The decoupling between different harmonics of X(a^, t) for small a is used for efficient 
computation of the second-order statistics of the time-periodic system (3.3). 



Appendix D. Computing the modifications k2 and 62 to k and e 

We next show that the averaged effect (over one period T) of fluctuations around the 



mean velocity on /c2 and 62 can be obtained from Xo,2(^)- Following (3.9), the second 



order correction (in a) to the auto-correlation operator of velocity fluctuations v averaged 
over one period T, (1/T) (v( • , a^, (g) v( • , a^, t)) dt, is given by C(a^)Xo,2('^)C'+(a^). 
The kinetic energy of fluctuations around the base flow and its rate of dissipation are 
given by (cf. equation ( |2.13 )) 



/ Kk{y,y,n)dn, e2{y) = / K^{y,y, n) dn, 

J K J K 



where Kk{y^£,^hi) and K^{y^^^K) are the kernel representation of the operators Nk and 
Xe, respectively 

Nk{K) = (l/2)(C,Xo,2C+ + aXo,2C+ + C^Xo,2C+), 

iK.x dy Cu Xo,2 C+ + Kxi^z Cu Xo,2 C+ + mz Cy Xo,2 C+ S+) + 
dyCuXo,2Cid^ + S^C^Xo,2C+5+ + /^^aXo,2C+ + 
hc^ Cyj Xo,2 + hil Cu Xo,2 . 
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Appendix E. Computing the effect of fluctuations on the mean 

velocity and the required power in the flow with control 

Here, we show how the effect of fluctuations on the turbulent mean velocity is deter- 
mined by the modification to the turbulent viscosity, ut2^ in the flow with control. The 
second-order correction to the mean streamwise velocity, [/2, is obtained by substituting 
U and Px from (3.12) into (2.4), and collecting the terms quadratic in a 

(^^T2(y)c/o(y) + 42(y)f^o(2/)), 



which yields the fohowing solution for U2 

U2{y) = -P.aMv) - 



1 



(El) 



The second-order correction to the driving pressure gradient, Px,2i is obtained from the 
requirement that the bulk flux of the flow with control remains constant, i.e. 



/: 



U2{y)<iy = 0. 



Integrating U2{y) in (El) from —1 to 1, and enforcing the above requirement yields 



'2UbJ-i J-1 1 



d^dy. 



An equation for Wp^2 is determined by substituting Vt from (3.11) into (2.7) 

{i + i^To{y))w;[2{y) + i^Toiy)w;,2iy) - mrcvtWp,2{y) = 
-My)w;\o{y) + 42(y)w^;o(y))- 

Finally, the fourth-order correction to the required power, %nreq,2, is obtained by sub- 
stituting Wp from ( |3.12[ ) into ( |2.1Q| ) 

_Ityi I IV' ^1 ^ 



%n.e„2 = ioo^im(i^;,|^_^ - w;^2\,^;). 



Appendix F. Computing the dominant characteristic eddies 

The most energetic flow structures are obtained from the principal eigenfunctions of 



the averaged (over one period T) auto-correlation operator, Xo(a^), defined in § 3.1 

[Xo{K)<p{-,K)]{y) = X{K)<p{y,K). 

At each k, A(k) is the largest eigenvalue of Xo(a^), and (f){ - ^ k,) = [v 7^]^ is the corre- 
sponding eigenfunction. The magnitude of (f) is determined by the requirement that the 
kinetic energy of fluctuations associated with cj) is equal to A(a^). On the other hand, the 
phase of <p is determined by requiring spatial compactness of v{x^ z) around x = z = 
in the lower half of the channel ( Moin fc Moser|p.989 ). This is achieved by making sure 



that J^^v{y^K) dy is a positive real number for all k. We note that enforcing compactness 
on u{x^ y^ z) yields similar results. 

In the uncontrolled flow, the principal eigenfunctions (f){ - ^k) account for approxi- 
mately 29% of the total kinetic energy; compare Kxi^z^if^) shown in figure fl6] with 



KxI^zE{k.) shown in figure 12(a) Furthermore, the two largest eigenvalues of Xo{k,) are 
almost equal to each other. As shown in figure \T6\ the difference between them is negli- 
gible for the values of k that correspond to the most energetic modes of the uncontrolled 
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Figure 16. (Color online) Premultiplied largest eigenvalue, KxKzX{k), of the auto-correla- 
tion operator Xo{k,) in the uncontrolled flow with Rr = 186. The contours show the rela- 
tive difference between the premultiplied two largest eigenvalues of Xo{k). The contour levels 
{10~^, 10~^, . . . , 10~^} decrease in the direction of the arrow. 



v{y, Kx — 2.5, Kz — 6.5) r]{y, Kx — 2.5, Hiz — 6.5) 




,41 , , . 1 -u.x^i ^ ^ ^ 1 

-1 -0.5 0.5 1 -1 -0.5 0.5 1 



(a? (bl 

Figure 17. (Color online) The real (black) and imaginary (red) parts of the principal eigen- 
functions {v^ (a); 77, (b)) corresponding to the two largest eigenvalues of the auto-correlation 
operator Xo{k) for the most energetic mode hix — 2.5 and hiz — 6.5 in the uncontrolled flow 
with Rr = 186. The two eigenfunctions are differentiated by solid and dotted curves. 



flow. In fact, the eigenfunctions corresponding to the second largest eigenvalue account 
for almost 28% of the total energy. This indicates that examining the effects of the eigen- 



function corresponding to the second largest eigenvalue is equally important. Figure 17 



illustrates that the eigenfunctions corresponding to the two largest eigenvalues of Xq at 
Kx = 2.5 and hc^ = 6.5 are equal to each other in one half of the channel and are mirror 
image of each other in the other half. Imposing the ^-compactness criterion on these two 
eigenfunctions aligns them in the lower half of the channel and it mirror-images them in 
the upper half. This implies that the flow structures that are obtained from the principal 
eigenfunction of Xo{k,) in the lower half of the channel account for approximately 57% 
of the total kinetic energy of fluctuations. 

The velocity components and w in the wavenumber space are obtained by acting 
with the operators C^, C^, and on the principal eigenfunction of the operator Xo{k,); 
see (3.4) for the definition of these operators. We use the flow symmetries in the span- 
wise direction ( Moin fc Moser|p.989 ) to determine the velocity profiles for the dominant 
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characteristic eddy in the physical space 

u{x, y,z) = 4 Re {u{y, k) e^'^^^) cos{k,zz) dn, 

v{x,y,z) = 4/ Re(v{y,K)e"^''^)cos{hizz)(^i^, 
= -4/ ImK,,.)e-)sin(...)d.. 
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